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P I ■ Abstract 

■ It has been found in numerical experiments [14] that when one removes a sector from an 

elastic sheet and glues the edges of the sector back together, the resulting configuration is 
radially symmetric and nearly conical. We make a rigorous analysis of this setting under 
two simplyfying assumptions: Firstly, we only consider radially symmetric configurations. 
Secondly, we consider the so-called von-Karman limit, where the size of the removed region as 
well as the deformations are small. We choose free boundary conditions for a sheet of infinite 
size. We show existence of minimizers of the suitably renormalized free energy functional. As 
a by-product, we obtain a lower bound for the elastic energy that has been conjectured in the 
related context of d-cones [16]. Moreover, we determine the shape of minimizers at infinity up 
(^■^ \ to exponentially decaying terms. 
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1 Introduction 



' The folding of paper or other thin elastic sheets is one of the many examples of energy focusing 

in the physical world. Starting in the late '90s, there has been a lot of interest in this problem 
in the physics community [14, 3, 5, 18, 6, 10, 13, 15, 4]. In particular the crumpling of paper 
(i.e. the crushing of a thin elastic sheet into a container whose diameter is smaller than the size 
' of the sheet) which results in complex folding patterns has drawn a lot of attention. It has been 

conjectured that the energy density per thickness h of such a folding pattern scales with h^/^. One 
major contribution in the rigorous analysis of this problem is [9], building on ideas from [17]. 
Here wc focus on approximately conical deformations of thin elastic sheets, that can be viewed 
as (one kind of) building blocks of crumpled deformations. Consider for example a sheet that 
is pushed into a hollow cylinder, such that the indentation of the sheet is small. The resulting 
structure is called a d-cone (developable cone). In the physics literature, it has been discussed 
e.g. in [3, 5, 14, 18]. There are several remarkable features of the d-cone, one of which is that the 
tip of the d-cone consists of a crescent-shaped ridge where curvature and elastic stress focus. In 
numerical simulations it was found that the radius of the crescent -Rcrcs. scales with the thickness 
of the sheet h and the radius of the container i?cont. as i?cios. ^ '^^^'^^cont • This dependence on the 
container radius of the shape of the region near the tip is not understood [18]. As argued in this 
latter reference, it cannot be explained by an analysis of the dominant contributions to the elastic 
energy, which are: The bending energy from the region far away from the center, which is well 
captured by modeling the d-cone as a developable surface there; and the bending and stretching 
energy part from a core region of size 0{h) where elastic strain is not negligible. The result of this 
(non-rigorous) argument is an energy scaling E ^ h^{Ci\\ogh\ + C2). This is a natural guess - 
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the situation here bears some resemblance to vortices in the Ginzburg-Landau model, where this 
is the right energy scaling [1]. 

In [16, 2] the scaling of the elastic energy of a d-cone with respect to its thickness h has been 
analyzed in a rigorous setting. The result from [16] is 



The lower bound does not achieve the conjectured scaling behaviour, and it seems that this claim 
can not be proved with the methods used in [16]. One of the aims of the present paper is to present 
a setting that is similar to the d-cone, but somewhat simpler from a technical point of view. We 
will show that a scaling result analogous to the one above without the log log h terms on the left 
hand side holds in this simpler setting. 

Crucial for a better understanding of the energy scaling in the two-dimensional model is some 

control over the rate at which the deformed configuration converges to the conical one when one 
goes in radial direction from the center of the disk towards its periphery. For general deformations 
of two dimensional sheets, it is very difficult to obtain this information. Here, we investigate a 
model for which it is natural to consider only radially symmetric configurations. This simplifica- 
tion allows us to quantify the decay of the difference between minimizers of the elastic energy and 
conical deformations for large radii. 

The idea is the following: To create an approximately conical deformation of an elastic sheet, 
we cut out a sector of angle /3 and glue the edges of this sector back together. We consider the 
minimization problem for the free elastic energy within the class of radially symmetric configura- 
tions - this seems to be a natural restriction of the set of allowed configurations. In particular, 
we will investigate the behavior of minimizers within this class as the distance r from the origin 
in the reference configuration tends to infinity. (We consider a sheet of infinite radius with free 
boundaries.) 

Apart from the restriction to radially symmetric configurations, wc make one more simplifica- 
tion in comparison to the situation in [16]: We use the so-called von-Karman approximation of 
non- linear elasticity [7, 11]. This means that the size of the removed sector /3 and the out-of-plane 
component of the deformation arc supposed to be of the order of some small parameter e, while 
the in- plane deformation is of order All terms in the elastic energy of order e'', k > 4, and of 
order h'^e'^, k > 2 are discarded. 

As we will explain in Section 2, these considerations lead to the definition of the free elastic 
energy density 



where \ = h/e and the deformation of the sheet is given as a map from spherical to cylindrical 
coordinates by 



/i2(C7i I log h\ - C2 log I log h\)<£h< h\Ci I log h\ + Cs) . 



(1) 





(2) 



with W' = w. The renormalized energy functional is 



W -J>IRU{-hoo} 



(3) 



{u, w) ^ lim^^oo Jo"" rdr (pf (r) - X^^^^) 



where 



W = {{u,w)G W/„'<f ([0, oo), : u{0) = w{0) = 



} 
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and ijj is some cutoff function with ip{r) =0 for r close to and ip{r) = 1 for r > 1. The boundary 
conditions included in the definition of W assure that the energy density is integrable in a neigh- 
borhood of the origin for the deformation eq. (2). 

The aim of the present contribution is to prove 

Theorem 1. The functional Ex from eq. (3) is well defined and bounded from below. It possesses 
minimizers in W (for which E\ takes on finite values). Furthermore, each minimizer («, li) satisfies 

u{r) =- — h o(exp(— (T-\/r/A)) 
w{r) =1 + o(exp(— cr-\/ r/A)) 

as r — )• oo for any a < 2y/2. 

As a side product of the proof of Theorem 1, we will get a lower bound for the elastic energy 

when the radius of the clastic sheet in the reference configuration is assumed to be finite. This 
lower bound is better than the analogous one from eq. (1) in that the log log-terms are not present. 
To give an idea how this "improved" lower bound comes about, let 

h= [ rdrpf . 
Jo 

The first step to establish the lower bound in the present setting is the right renormalization of 
the clastic energy density. We expect a logarithmic divergence in A of X~^I\ as A ^ 0. Thus we 
make the replacement 

The key step is now to find a change of variables that makes it obvious that 

/,*(,..M-A«) (4, 

is bounded from below by some constant times A^. As we will see in Section 3, such a change of 
variables does exist, and will leave us only with manifestly positive terms in the renormalized energy 
eq. (4) plus some divergence- like term that will be estimated in a suitable manner in Proposition 
1. Thus we get the sought-for lower bound 

A-'/A>|logA|-C. (5) 

This paper is organised as follows: In Section 2, we motivate and define our model. In Section 3 
we prove the existence of minimizers of the elastic free energy functional. The lower bound for the 
elastic energy will be a corollary of one of the preparatory propositions. In a remark at the end of 
that section, we will discuss a pathology of the model presented here. In section 4, we show that 
minimizers converge to the conical configuration at infinity. 

2 The model 

Cutting out a sector, glueing the edges back together. For small e > let be defined 

by 27r/(27r - /3(^)) = y/t+e^, and let 

= |R2 \ |(xi,X2) : X2 < < a;i, -/^^^^ < arctana;2/a;i < o| 

We define a deformation y : B^^^ — )• whose image has rotational symmetry. We are going to use 
cylindrical coordinates {r,ip), 

y{r,ip) = U{r)eif^ +V{r)e, (6) 
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where 



U : [0, oo) [0, oo) 

V : [0,oo) [0,oo) 

e(^) = (cos(^(^),sin(p(^),0) 



We calculate 



Vy^Vy - Id = ((C/')^ + (V)^ - l) + (^^2^^^' - 1^ 

V^y = (l7"e(^) + y"e^) O O + VTTe^ e^f) O (e<^ (g) + (g) e^) 



Definition of the elastic energy. As a starting point, we choose the elastic energy density to 
be 

Pel. =|Vy^Vy-Id|' + /i2|vV 

where h is a parameter for the thickness of the sheet under consideration. The first term on the 
right hand side represents the stress energy density, the second one the bending energy density. 
This is a standard ansatz for the elastic energy, for a justification see e.g. [9]. 

Von-Karman ansatz, change of variables. Now we mal?;e the ansatz 

U{r) =r + 0(e2) 

V{r) =0{e) , (7) 
where e is some small parameter. The following changes of variables will be convenient, 

U{r) =r + — {u{r) — r) 

V'{r) =sw{r) (8) 



and, alternatively. 



Uir)=r+^u{r) + X^^-r 



V'ir) = e{iP{r/X) + w{r)) , 

where A = h/e, tp € C°°([0, oo) with ^p{r) = near r = and ip{r) = 1 for r > 1. 
A deformation of finite elastic energy must satisfy the boundary conditions 

u{0) = u(0) = w{0) = w{0) = . 

A short computation shows that the elastic energy density is given by 

- \ 2 / ^2 

~2 



6\ I /n('^6\2^ 
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We build our model only considering the leading terms in e for fixed A, defining 



/'A =li™;£ Pel. 

= i{<b'-l + u'f+(l)\x'U'' + ^]\. (9) 



Renormalization and rescaling. The (u, w) variables have been chosen such that we expect 
u, u', w, w' to vanish as r — )• oo. By an inspection of the energy density Px'{r) we expect that the 
integral 



/ rdrpf(r) 
Jo 



diverges logarithmically as R ^ oo. To have sonic hope of a meaningful limit for R — )• oo, it seems 
a good guess to introduce the renormalized functional 



{u,w)^ j rdr (^pf - {r) - \^ 



(10) 



This definition is to be understood in the sense that E^{u,w) = oo if the (Lebesgue-) integral on 
the right hand side does not exist. In the sequel, we will set A = 1, and derive all results for this 
value of A. The general case can be recovered by the change of variables r r/X, which we will 
do at the very end. In fact the change of variable formula yields 

E^{u, w) = X^E^'^ {X-^u{X-),w{X-)) (11) 

for any u, w G Wj^'^. 

We will use the following notation: 

Ef =E^ 

E^{u,w) =E^{u,w) (12) 

3 Existence of minimizers 

Proposition 1. For the functional E^ : W ^ IR defined in eq. (12), 

E^{u,w) > -C, 

where C is some numerical constant that is independent of u,w and R. 

Proof. We bound from below the integral from [0,1) and [1, i?] scperately. The contribution from 
[0,1) is trivially bounded from below since p^^- = pf- is positive. We set 



Xi = Xi(u, w)= I rdrp"^- > . 



Obviously, 

E\u,w)>Xi-C , (13) 

where C = tp'^dr/r < oo. 

Next we consider the contribution from [1, .R]. Note that for r > 1 the (renormalized) elastic energy 
simplifies to 

p^'-(r) - ^ = (2.0 + + nO V + - i (^)' + _L . (14) 
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To show boundedness below of the integral /j^rdr(p^^' — l/r^), we will show that the absolute 

value of the integral dr{u/ry can be bounded by square roots of the manifestly positive terms 
in -u;). Obviously, 



u{R) 
R 



-u{l). (15) 



First we note 

nl 



2 



2 / driiu' 



< 2 ( / dr — 1 ( / rdr-'2 



1/2 / rl \ 1/2 

u" 



Jo r , 

(pi -2\l/2 / / /-l /-l W"'^^^ 

J dryj i^i^J rdr {w^ - 1 + u'f + J rdr{w^-l fjj (16) 



In a similar way, we estimate the supremum of w on [0, 1], 



Furthermore, 



sup w{xy<2{ dr — / rdrw'^] <2Xi (17) 

a:e[o,i] \Jo r J \Jo ) 

[ rdr (w^ _ 1)2 < max f 1, sup ui{x)A < AX^ + 1 . (18) 
Jo \ xe[o,i] J 

Inserting eq. (18) into eq. (16), we get 

/ /■! f,2X 1/2 / rl N 1/2 

|u(l)|2 < 2\/2(y dryj i^J rdr{w'^ -1 + u'f +4X^ + 1 j 

and hence 

Ki)|=Ki)-i/2| 

<C{X^^^ + 1) (19) 

We turn to the estimate of u{R)/R. We may assume u{R) > |w(l)|. Choose a such that u{a) = 
|u(l)|, u{r) > u{a) for r > a. Let furthermore 



A = < r e [a,R] : u{r) = max u{s) 

1_ sS[a,r] 

Ai ={r e A : u'{r) < 2} 
Aj ={r e A : j < u'{r) < j + 1} 
A200 = Uj>2 Aj 



We have u{R) — u{a) — u{R) — \u{l)\ > 0, and the following case distinction: Either 

/ u'{r)dr > (case (i)) 

JAi 2 



or 



f . , . u{R)-\u{l)\ 

/ u (r)dr > — — (case (u)) . 



In case (i). 
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and 



dr—>- dr(|u(l)|+2r) 



1 

R 

C 



-r^ + 2r'^\u{l)\+r\u{l)\^ 



iu{R)-\u{l)\)/4 




>^(Ki?)r-Wi)H , 

where we have used Young's inequality to estimate the mixed terms in u{R), and C is some 

numerical constant. Thus we get in this case 



where 



In case (ii), 



X3 = X3{u,R)= / —dr. 

Ji r 



u{R)-\uil)\ 



2 <^(^ + l)|A,| 

i>2 

<3^(z-l)|A, 



i>2 



3 / EXA^(^-I) 

•'A ,.^o 



i>2 

< 3 / drlw^ -l + u'\ 
/a 



< 3 ( y rdr|w;^ -1+u 



J dr/i 



< C 



A 

R y/2 

rdr (2w; + + u'f (log i?)^/^ 



Thus, choosing R such that i?^/^ > log-R, we get for case (ii), 

R-^^\{R) < C(xf ^ + X^^^ + 1) 

where 

X2=X2{u,w,R) = J rdr {2w + + u'f . 

By eqs. (20) and (21), 

r-'/MR) <c[i + xf/^ + xy + xy^) 

We put together eqs. (13), (14), (15), (19) and (22) to obtain 



/' 

JO 



rdr 



>Xy + X2 + Xs 

- c (xf * + xy^ + xij^ + xy^ + 1) 



(20) 



(21) 



(22) 



(23) 



This proves the proposition, since the right hand side is bounded from below as a function of the 
Xj, i=l,2,3. □ 



Corollary 1. There exist constants C^^C^ such that 

I log A| - Ci < inf 7a < | logA| + C2 . 

Proof. By eq. (11), 



Using eq. (12), we get 



El{u,w) = X'E^ (Au(-/A),^-/A)). 



iuiEl < X^E^ '(0,V) < CA^ 



(24) 



which proves the upper bound by the definition of Ix (we recall I\ = E\ + }? J^{'il){r/X)'^/r)dr). 
The lower bound follows from eq. (24) and Proposition 1. □ 



We next send i? — > 00 in the definition of E^. For {u, w) G W, we define, 

, 2 



E^{u, 111)= rdr 



+ / rdr 



(25) 



This definition is to be understood in the sense that E~^{u,w) = +00 if the integral on the right 
hand side of eq. (25) does not exist. Note that the derivation of eq. (23) also implies 



/' 

Jo 



rdr 



<Xi +X2 + X3 

+ c (xl'' + xl'^ + xl'^ + xl'^ + 1) 



(26) 



and hence Mmn^oo rdr(p'^'' — V'^/?'^) = 00 if and only if E~^{u, w) = 00. 

For {u, w) gW such that either the limit limr_>cx> exists or E~^{u, w) = 00, we set 



E{u, w) 



Hfcl + uil) + i - rdr|^ if lim^^oo ^ exists 

if E~^{u,w) = 00 

We will see in Corollary 2 that E is in fact well defined on W. By eq. (14), 



E^{u, w) — \im.r 
00 



E{u,w) = lim E"-{u,w) 



(27) 



(28) 



if lim„ 



u{r) 



exists. 



Lemma 1. Let k > and g : [0, 00) — )• IR 6e continuous with \g\ > k. Then there exists a constant 
C{k) > such that 



lim inf 



1 



inf 



Proof. By the change of variables x = yl, f{y) = f{yl), the claim is equivalent to saying that the 
infimum of 



is bounded below by some positive number as ^ — )■ 00. Assume this were not the case. Then there 
exist sequences {lj}j C R, {fj}j C M^^'^([0, 1]) such that 



I' ^{fj+ijaY + f^dy^o 



(29) 
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Hence / /J — 0, and for every 6 > there exists jo such that for j > jo there exist yj^i € [0, 1/4], 
yj,2 e [3/4, 1] such that \f{yj,i)\ < 5 for i = 1, 2. Thus 



yj.i 



and 



This imphes 



/ (fjiy) + lj9]dy 



< 26 



> -Lk-25. 



hminf (/'(y) + IjgVdy > Cn^ 



for some numerical constant C > 0, which is a contradiction to eq. (29). 
Lemma 2. Let {u,w) G W and E~^{u,w) < oo. Then 

lim w(r) =0 

r— >oo 



-1/3 



|ii(r)| <C 



□ 



Proof. Note that only depends on and w'"^ which implies E^(u, w) = E^{u, \w\). Hence we 
may assume w = w + V' > 0. Let w{t) = w{e'-), u{t) = u{e^), 5 > Q and assume 



lim sup w {t) > 6 



(30) 



We have liminft_,.oo w'^{t) = since otherwise | 2w+w'^\ oo. Thus, we can pick two sequences 
{sj}j£iN, {tj}je^ C [0,oo) with the following properties: 



Sn < tn < s 



Sj , tj 



w^{sj) = 6 
^Hh) =6/2 

6/2 < w{t) <6 for t e [sj,tj] 

Note that by the last inequality, we have 

|e* {2w{t) + w^it)) I > C{6) for t e [sj,tj] . 



Now 



and hence 



6/2 < J^' {w^ydt < 2 w^dtj w'^dt^ 



f 



w^dt > 



f*^ w''^dt 

J Si 



^ 00 as J ^ 00, 



(31) 



since the denominator goes to by the integrability of w'"^. As a consequence, \tj — Sj\ ^ oo. 
Thus (by a change of the domain of integration) the conditions of lemma 1 are fulfilled for f = u, 
g{t) = e*(2«}(i) + w)^(i)), k = C{6), and we get 



p {e\2w + w'') + u'Y +u^dt 

J Si 
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which is a contradiction to E'^{u, w) < oo. This proves the first statement. 

Let M = supw^, and ti,t2 be such that w{tif = M,w{t2f = M/2, M/2 < 'up-{t) < M for 
t € [ti,f2], where for definiteness we have assumed ti <t2- Now we estimate 



M/2< f\wydt 
Jti 



1/2 ^ 1/2 

<2( / ^i^di] ( / w'Mt 



tl 

t2 \ 1/2 



With obvious changes, the estimates also hold true for ti > ^2- Assuming M > 2, we may apply 
Lemma 1 with f = u, g = e*{2'w + w-^) and k = 1. This yields \t2 — ti\ < C E^{u,w), and the 
second claim follows. 

The third statement only has to be proved for negative u{R), for positive u{R) it has already been 
shown to be true in the proof of Proposition 1. 
Choose to such that 

f-OC POO 

/ dtu'it) = / dr(u2(r)/r) < 1 

Then for i > to, J.~^^ u'^dt < 1 and hence min(g[j^j_|_i] \u{t)\ < 1. We may then choose 

ti e argmint£[i^i+i]|u(t)| , 

and set := c*'. For given R, choose i such that ri < R < rj+i. Now we proceed in almost the 
same way as in the proof of the proposition. Let 

A = < r e [R, rj+i] : uir) = max u(s) > 
I selR,r] J 

Ai ={r e A : u'{r) < 2} 

Aj ={r e A : i < u'{r) < j + 1} , j>2 

A200 = Uj>2 Aj 

We may assume — 1 > 0. The case distinction is this time 

/ u'{r)dr > '"^-^^l ~ ^ (case (i)) , 
u'{r)dr > ~ ^ (case(ii)) 



/ 

J A, 



In case (i), 

|A,| > 

4 

and hence, by a calculation just as in the proof of Proposition 1, 

rn+i „2 



dry>C - 1) . 



IR 

In case (ii), 

(R \ ^1"^ 

rdr{2w + w'+u')A 

Hence \u{R)\R~^/^ is bounded, and the lemma is proved. □ 
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Corollary 2. The junctionals (3) (for \ = 1) and (27) are well defined and hounded from below. 

Furthermore, E{u,w) ~ E{u,w) for all {u,w) G W. 

Proof. Below eq. (26), we have already remarked that limii_,.oo J^rdr{p'^^- — ip'^/r'^) = oo if and 
only if E~^{u,w) = oo. By Lemma 2, E~^{u,w) < oo implies liuir^oo u{r) /r 0. The claims of 
the corollary then follow from eqs. (27), (28) and Proposition 1. □ 

Proposition 2. Let {{uj,Wj)}j^t^ CW be a minimizing sequence for E, i.e., 



lim E{uj,Wj) = mi E . 



Then 



(i) there exists a subsequence that converges to some (uo,wq) in Wi^'^((0, oo)). 
(a) lim^^o ■wo(r) = limr^o w'o(?') = 
(Hi) E{uo,wq) =va.iY^ E. 

Proof. (!) By Corollary 2, E is bounded from below, so infw E > — oo. 

Let < a < 6 < oo. We will show that E is coercive on VF^'^([a, 6]). We have E{uj,Wj) < 
C for some constant C by the minimizing property of the sequence {{uj,Wj)}j. This implies 
E^{uj,Wj) < C. Hence the boundedness 



/ u]+wfdr <C{a,b) 

J a 



(32) 



is clear. In the sequel, we drop the subscript j. Now 



/ u'^dr <2 / dr(u;2 + 2wV + V'^ - 1 + u'f + {w^ + 2wV + V'^ - 1)^ 

J a J a 

<2 ^a-^£;+(u,w) + j dr{v? +2wi) + i)^ -if^ 



(33) 



We will estimate the second term splitting the integral into two parts, with integration domains 
[a, 1] and [1,6]. Firstly, using eq. (18), 



/ {v? + 2u;V + V'^ - 1)^ < max (sup |^^;|^ l) < A{E+{u, wf + 1) . 

J a 



(34) 



For r >1, up' + 2w^ + — 1 = 2t« + uP. We may assume w{r) > — 1 (cf. the proof of Lemma 2). 
Let w{t) = w{e*),u{t) = u(e*) and 



2w{s) + w{sfds 



Now we are going to use the interpolation inequality 

= ([0,oo)) 



||2?D + w^|||2( 



W'\\hi[0,oo)) < C\\W\\L.i[0,oo))\\W''h^[0,oo)) 



{35) 



We will estimate the terms on the right hand side. By lemma 2, sup < C {E{u,wy^'^ + 1). 
Thus 

\\W"\\l2 = \\{2 + w)w'\\l2 <C{E{u,wy/'^ + 1)\\w'\\l2 

<C{E{u,w) + 1) 



Also, 



\W{T)\ < 



f 

J OO 











+ 








J OO 
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Using Holder's inequality the first term on the right hand side can be estimated by 

We integrate the second term by parts and get 

lim Te-^MfT) - e-'"M(r)l + / dte-%(t) 
By Lemma 2, limi-^oo e~'^u(r) = 0, and again by Holder 



J oo \J oo J 



1/2 



Putting the last four equations together, 

\W{T)\ < Ce-^ +m(T)) 

and thus 

\\W\\l2 < CEiu.wf/'^ 
Inserting the estimates of ||T^||j;,2, ||T^"||i,2 in the interpolation inequality (35) yields 

/oo 
{2w + w'^fdr = \\2w + w'||i2([o,oo)) < C{E+{u, w) + l). 

Inserting this result and eq. (34) into eq. (33), we get 

< C(a) {E+{u,w) + 1) . 

Since 2 + u;>2 — ■^>1, eq. (36) implies 

/oo 
w^dr=\\w\\l,^^^^^))<C{E+{u,w) + l). 



Also, 



and thus 



J a 



W 



< sup < C{E+{u, w) + l), 



l-b 

/ w'^ <C {E+{u,w) + l) . 

J a 



(36) 



(37) 



(38) 



(39) 



Reintroducing the subscripts, eqs. (32), (37) and (39) show the boundedness of {{uj,'Wj)}j in 
Wi'^([a, &]) for fixed a,b. We now pick a subsequence {iu^^\wj^^)}j that weakly converges in 
^^'^([1/2, 2]), and inductively we choose subsequences {iuj''~^^\wj''~^^^)}j of {{uj''\wj''^)}j that 
weakly converge in W^''^{[l/{k+l),k + l\). Let (uq, wq) be the weak limit of the diagonal sequence 

{{uf\wf)}j. 

(ii) Let K > be such that ?/) = on [0, k]. Then u = u,iv = w on [0,k]. We introduce the normed 
space 

Wi := ju; e W^^^{0,k) : ||w||i,2((o,„),dr-/r)) + lk'IU2((o,K),rdr)) < ooj 
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with norm ||w|lwi = l|u'llL2((o,K),dr/r)) + l|w''llL2((o,K),rdr))- {{w,w') : w G Wi} is a closed subspace 
of L^(dr/r) x L^(rdr) and hence Wi is reflexive and bounded sets are weakly sequentially compact. 
For the diagonal sequence from (i), we have 



I 



rdr iwf + < E+{u, w) 



and hence wj wq in Wi. In particular, wq € W^. Thus 

rdrw'^ + ^-)-Oas^^O, 



Jo 



and 



\ 1/2 / ^ ^\ ifi 

^ rdrw'A n rdr^ 



as ^ ^ 



which proves the claim lim^^o '^o('') = 0. Analogously one shows limr_^o Uo{r) = 0; the bounded- 
ness of rdru'^ in this case follows from the estimate 

/ rdruf <2 / rdr {u'j + «;| - l) + (wf - 1)^ 
Jo Jo 

<2E+{uj,Wj) + K^(sup \wj\'^ + 1) 
<C{l + E+{uj,Wj)+E+{uj,Wjf) . 

(iii) We have Uj — ^ UQ,Wj — ^ Wq in Wj^'^(0,oo), and hence {uj,Wj) — >■ {uq,Wo) uniformly on 
compact subsets of the positive real line. Hence, by Fatou's lemma, 

lim'mf E~^{uj,Wj) > E~^{uq,'Wo)- 

j->oo 

Also, limj^oo = uo(l)- Hence 

E{uo,'Wq) = E~^{uq,wo) +uo{l) — I < liminf £'(u,-, Wj) 

Jo r j^oo 

which proves the claim of lower semicontinuity. □ 

Remark 1 (Self-penetration of solutions) . The von Karman model displays a pathology at the origin 
for the situation we want to model. Namely, the solutions we have found show interpenetration of 
matter. Condsider again the Euler-Lagrange equation obtained by variation of u, 



(r {w^ -1 + u'))' =-. 



u 
r 

Since w — )• for r 0, the qualitative behaviour of solutions u near the origin is the same as the 
one of solutions of the linear equation 



(r(«'-l))' = ^. 

r 

The solutions of this latter equation are given by 

^rlogr + Cir + C2r~^ 

The integration constant C2 has to be set to zero to fulfill the boundary condition ?i(0) = 0. Going 
back to eq. (8), we see that the value of U will be negative in some punctured neighbourhood of 
the origin and we have sclf-pcnctration of the solution (somewhere in the region r ^ exp (— e^^)). 
We expect that this pathology could be cured by including nonlinear or higher order terms in u 
in our model. We refrain from doing so, since the main aspect of this work is the analysis of the 
solutions away from the origin. 
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4 Decay properties 

We now turn our interest to the decay properties of minimizers. First we show that the decay of 
the energy density is fast in a weighted L^-sense. This will level the field for an application of 
stable manifold theory, by which we prove that u and w decay exponentially. 

Proposition 3. Let {u,w) be a minimizer of the functional (27), and < a < 2. Then 

r.i+«dr (^{2w + + u'f + ^ + w'^^ < oo . 

Proof. 1. The idea is to integrate the Euler-Lagrange equations against suitable test functions. 
We may assume a > 1, since the claim is stronger the larger a is. Let 

• f 2-a 
A = min a 



Let ro > 2 be chosen such that 



25 

Ke(l/2,1). 



1 2?/! -I- up' 1 

Kr)|<-, ^(T^<iforr>ro-1 (40) 

which is possible by Lemma 2. Let R > ro. Now let ^ e C°°([0, oo)) with ^(r) = for r < ro — 1 
and ^(r) = 1 for r > ro, and 



Define the test functions 



5(r)r« if r < i? 



(j)[r) =M(r)xa,7?,(r) 

2w{r) + w'^{r) 
= 2il + w{r)) ^"'^^^^ • 

We claim that E^{u + t(p,w + tip) is finite for any t G R. To see this, we look at each of the 
summands of the integrand of the energy functional soperately for r > R, and estimate them from 
above by functions that are integrable with respect to the measure rdr: 

{2{w + t<p) + {w + tipfiu + t(t>)'f = {{2w + + u'){l + tXa,R) + iV^ + tux'^^n)^ 



<C{R,t,K) (j^2w + w'^ + u'Y + (2w + 



2\4 4 , " 



Since w is bounded by Lemma 2 and by the decay of Xa,n. at oo, {2w+w'^)'^Xa r i'^ ^^([Oj oo), rdr). 
Also, by Proposition 2, 2w + + u' and u/r are in L^{[0, oo), rdr) which proves that the left hand 
side of the above equation is in ^^([O, oo), rdr). Similarly, 

{u+t(py u'ji+txf ^ 

— ;:2 — = ;:2 ^ <^1^'*)72 



{w + t^f = [w' [l + 



tXa.R A , 1 \\ 2w + w2 



2 V (1 + HV/ 2(l + w;) 
<C{R,t){w''^ + 



,2 , {2w + w^f 
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where each of the three estimates imphes that the left hand side is in i^"'^([0, oo), rdr) by the results 
from Proposition 2. This proves the claim of the finiteness of E^(u + t(p,w + t(p), and hence of 
E{u + t(l),w + tip). The minimality of E{u, w) implies 



=—E{u + t(t),w + t(p) 



r°° ( f 2w + w'^ 

--C{ro) + 2j rdr I {2w + + u') \^2{1 + w) ^^^^^^ x«,fi + u'Xa,R + ux'a,R 



,, , w'(2w + w'^) \ 2w + w^ , ^ 



2{i + wy 



=C(ro) + 2 1 rdr {2w + + u'f + — + u;'^ 

Jrn . 



+ 



,,2 



2(1 + w) 

2 
2 



Xoi,R 



■, , w'^{2w + w^) 2w + w^ . 

{2w + W +U)UX^^^- .2 Xa.fl+ „,, ^ . Xc.fl 



2(l + w)2 ^"'^ 2(1 + w) 
Note that u + t(j) = u and w + tip = w on [0, ro — 1], whence 



(41) 



rdr 



ro — 1 



(2m; + + u'Y + ^+W 



Xa,R 



+ 



{2w + w^+ u')ux'r 



w''^{2w + w"^) 
2(1 + ^)2 



Xa,R + 



2w + uf 
2(1 + w) 



The only thing that will matter for us is that by the definition of Xoi,R, C'(ro) does not depend on 
R. We set 



/•OO 

Xi = Xi (ro, R)= rdr (2u; + u;^ + u'fxoc,R 

J To 

Xi = Xi{rQ,K) = / rdr —Xa,R 

Jro ^ 
poo 

X3 = Xsiro, R)= rdr w'^Xc.,r ■ 



The Xi, i = 1,2,3 are manifestly positive. We are going to show that the absolute value of the 
other integrals in eq. (41) can be bounded by terms that are at most linear in the Xi and terms 
that arc independent of R. 
2. By Young's inequality. 



^ rdr(2«; + w'^ + u')ux'„^r <| (^J rdr (^i2w + w'^+ u'f + Xa.fl 



where we used 



<-(Xi+X2), 



C'(r-)r" + ^(r)ar"-i if r < i? 
-Ki?«+'^r-(i+'") if r > 



(42) 



(43) 



which implies jx^ /{| < max(K, a)xa,R/r =< aXa,R/r for r > rg. Secondly, by our choice of ro. 



/•°° ^ u;'2 2«; + w;2 1 |- 

/ ^dr 7^Xa,fi <7 / rdrw Xa,R 

Jro 2(1+m;)2 ' a J^^ 



ro 



(44) 
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3. Finally, we have to estimate the integral 



2(1+.;) 



Splitting the integral into two parts, with integration domains (0,ri) and (ri,oo) with ri > ro to 
be chosen later, we can estimate its absolute value from above by 



5 f°° 

C{ro,ri) + -l r&r \w' {2w + w^)x. 



2^,,/ I 

a,R\ 



(45) 



(The factor 5/4 stems from eq. (40); also, we assumed n < R.) By eq. (43) and Young's inequality, 
the integral can be estimated by 



5 5 
- J rdr \w' {2w + w'^)x'a,R\ rdr 



1 25 
< -X3 H 

- 4 ^ ^ 4 



,2w + 10^ 

W Xa,R 



rdr 



/ 2w + 



Xa,R ■ 



(46) 



The second summand on the right hand side above will be estimated with an interpolation inequal- 
ity in a similar way as in the proof of Proposition 2. We set 

W{T) = f{w\t) + 2w{t))x^'\t)&t 

where w{t) = w(e*), x(*) = Xa,ii(e*)- Note that the definition of W depends on R and 

, „ 9 \ 2 
2w + w 



f°° ( ' 



Xa,R = \\W \\l2(^[Ti^oc:)) ' 



(47) 



where we have set e-^^ = ri. We also set u{t) = u(e*). For every 5 > (to be chosen later) there 
exists a constant C{d) that only depends on S such that 



\W'h.^[T,,oc)) < 5||W^"IU=([T„oo)) +C(5)||W^|U.([T„cx,)) 



(48) 



For a proof of this interpolation inequality, see e.g. [12]. Using the triangle inequality, integration 
by parts and the Holder inequality. 



\W{T)\ < 



< 



(X 



di(w2 + 2'ii) + e"'{t')x^/^ 
1/2 



-T 



dte-'u'x'^^ 



dte 



-2t 



1/2 



dtie\2w + w')+u'yx] + u{T)e-^x^l^{T) 



+ 



f dtuU-'x"^)' 

J CO 



(49) 



Assuming e* > ro, we calculate 



(e-*X^/^(i))' 



(f - l)e-*x^/^ ife*<i? 
-(l + f)e-*x^/2 ife*>i? 



and thus |(e-*x^^^ (*))'! < 2e-*x^^^- Inserting this into eq. (49) yields for T > Ti: 



\W{T)\<^e-^ \ dt{e\2 



1/2 



w + w^) + u'Yx] +2 / dtu^x 



Ti 



l/2^ 



+ 



u{T)e-^x'/\T) 
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and hence 



\W{T)\'^ < 2e 



-2T 



Xdt{e\2w + + u'f +u^]+ u\T)x{T) 



Ti 



Xdt{c\2w + w^) + u'f + 



With these preparations, we can estimate the L^-norm of W, 

/ W{T fdT <2c"2^i 

<2e-2^^ (Xi + X2) 
We turn to the estimate of ||W^"||i^2. For e* > ro, 

= 2(1 + w)w'x^'^ + {2w + w^){t'^)' 
<4 [w'"^ + {2w + 'iJp'f^ X ■ 

Hence 

p dt\w"\'' < 4.(^X3 + p W'^dt^ 

We insert eqs. (50) and (51) into eq. (48) and get 

WW'Wl. < 46 (^3 + llW^'ili^) + 2C{d)e-^^' (X, + X2) 
Now we choose 6 such that 



4(5 



1 

1 - 4(5 ^ 25 



and Ti such that e^^ > ro and 



2C{S) 
1-46^ 



-2Ti 



< A. 



Note that Ti and hence ri only depends on ro and a, but not on R. Eq. (52) becomes 

l|W^'lli=([T„oo)) < A(Xi + X2) + . 

Putting together this latter estimate with eqs. (45), (46), (47) and (53), we get 



"•"^ , w'(2w + w^) 

rdr Xa,_R 

2(1 + u;) 



(50) 



(51) 



(52) 



(53) 



< C(ro,a) + ix3(ro,i?) + ^(A(Xi(ro,i?) + X2(ro,i?))) (54) 



4. Inserting eqs. (42), (44) and (54) into the integrated Euler-Lagrange equation (41), we get the 
estimate 

>Xi +X2+X3- |^|(Xi + X2) + ^Xs + C{ro, a) + ^^3 + ^A(Xi + ^2)^ 

>^(Xi+X2) + Jx3-C(ro,a). 

This inequality holds for every R, which implies that the Xi{ro,R), i <E {1,2,3} are uniformly 
bounded in R. By definition they are integrals of functions that are monotone in R, so by the 
monotone convergence theorem the limit function is integrable and the value of the integral is 
identical to the finite limit lim/{_>.oo ^ii^o, R)- CH 
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Wc already know that lim^^oo w{r) = for any minimizer (u, w). Now wc show that minimizcrs 
actuaUy have exponential decay at infinity. We will use the following standard tool from stable 
manifold theory: 

Theorem 2 ([8]). Let sq e R, ^ e R"^" a matrix with k < n eigenvalues with negative real part 
and n — k eigenvalues with positive real part , F : R" x [s, oo) R" with the property that for every 
K> 0, there exist Sq € [so,oo) and <5o > such that 



whenever \x — x\ < Sq and s > Sq. Then there exist 5 > 0, s > Sq such that for \p\ < S and s > Sq, 
there exists a k- dimensional submanifold M[s) of R" containing the origin such that the initial 
value problem 



has a solution x : [s, oo) — )■ R" for any p e M(s), and the property 

\x{s)\ = o(exp(— (T.s)) as s ^ CO 

for any a > such that the absolute values of the real parts of the eigenvalues of A are all bigger 
than a. Furthermore there exists r] > independent of s such that, if p ^ M{s), then 



\F{x, s) — F{x, s)\ < k\x — X 



—x{s) = Ax{s) + F{x, s) , x{s) = p 



(55) 



|a;|L~(so,oo) >v- 



Proposition 4. For any minimizer {u,w) of the functional (27), 

\{u{r),u'{r),w{r),w'{r))\ = o (exp(— cr\/r)) 



as r 00 for any a < 2\/2. 

Proof. The Euler-Lagrange equations in the {u, w) variables for r > 1 read 



2(1 + w)r{2w + w'^ + u') ={rw')' 
{r{2w + w'^ +u'))' =-. 



(56) 
(57) 



We choose i?o > 1 large enough so that 



- < 1 + w{r) < - for r > i?o , 
which is possible by Lemma 2. In this region we may divide eq. (56) by 2(1 + w), 



(58) 



r{2w + w'^ + u') = 



1 




2(1 + u;) 



Here and in the sequel, it is tacitly understood that r > Rq. Eq. (57) becomes 




Inserting this into the first equation again, we get a fourth order equation for w only, 




(59) 



We make a change of variables, 



s = \fr w{y/r) 



w{r) 
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such that eq. (59) becomes 

{1 + w) 
This can be rewritten as 



2w + 



16s I V2s(l + «!) 



{sw'Y 
8s3 



w'-^\s) = -64u;(s) + g{x{s),.s) + h{x{s),s) 



(60) 



where x{s) = {w''^\s),w''{s),w'{s),w{s)), and £f : R'* x IR+ ^ IR contains the nonUnear terms in x, 
/i : X R+ R the linear ones with coefficients 0{s~^). More precisely, 



\-\-w_\ S S'^ i- + w\ s 

- 96w^ - 32w^ 



h{x, s)=- -v/^'^ + ^w^^' + ^v/ 



(61) 
(62) 



S S'^ 

In particular, for /:=<? + /i, and given k > 0, there exist S,5 > such that 

\f{x,s) - f{x,s)\ <k\x-x\ 

whenever \x — x\ < 6 and s > S. Additionally, we have /(O, s) = 0. Now we may rewrite eq. (60) 
as a system of first order equations. 



ds 



x{s) = Ax{s) + F{x, s) 



where i^^ : R"* x R H> R" is given by F{x, s) = {f{x, s), 0, 0, 0) and 



-64 

Id3x3 



The eigenvalues of A are 4(±1 ±i)/-v/2. I.e., A has two eigenvalues with positive real part and two 
with negative real part. We are now going to show that for a minimizer {u, w) of the functional 
(27), the function w[s) = w(s^) satisfies 



/•OO 

/ ds|u;(4)(s)| 



< 00 , 



(63) 



where 5o = ^/Rq. Prom eq. (63) it follows that lims_,.oo (s) exists and is finite. By Lemma 2, 
«; ^ as s ^ DO. Hence 

hm y/^\s) = lim w"(s) = lim w'{s) = 0. 

S— >-00 S— )-0O S— 7-00 

It follows that a; — >■ 0. From Theorem 2, it follows that there exists s such that x{s) G M{s), and 
hence \x{s)\ = o(exp(— as)) for a < 2\/2. 
It remains to prove eq. (63). 
By Proposition 3, 



r r^/^dr{2w + w^+ u'f = 2 T s^ds (s{2w + w^) + =^ 

where u(s) = u(s^). In the new parametrization, eq. (56) reads 

1 



< DO , 



(64) 



2(l + u;)s^ si2w + vt) + ^ = -{sw') 



19 



The left hand side is square integrable by eq. (64) (the factor of 2(1 + «;) is irrelevant by eq. (58)). 
Thus 

r{sw"fds <C ( risu/y^+n/^ds) 

J So \JSo J 

<C I / r^/^dr{2w + w'^ + u'f + / r^'^drw 

\J Rq J Rq 

/ POO POO \ 

<C / r5/2dr(2«; + w2 + u')2 + / rdrw'M 

\Jro Jro J 



<00 . 



where in the second inequality, we have used the trivial estimate r > r^l"^ for r > Ro- This implies 



poo 

/ \m."\ 

J So 



1/2 



So 



1/2 

{sw")'' ] < 00 



(65) 



and the existence and finitcncss of the limit lims_>.oo w'{s) follows. Moreover, since by Lemma 2, 

lims^.oo w{s) = 0, we also have 



lim w^{s) = and sup |w'(s)| < C . 

s>So 



(66) 



In the new parametrization, eq. (57) reads 

^ ' „(3) 



s 8{l + w) 
which implies (using eqn. (58)) 



w 
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s 1 + w s(l + w) J ' 



< c 



u{s) 



l + w 



s{l + w) 



(67) 



Again using eqns. (58) and (66), we get 

rSi 



Jso 





u{s) 


\-fso 




s 



+ |u;"P + |w;'|2 



(68) 



for any 5*1 > 5*0, where C does not depend on Si. By interpolation, for any S > there exists a 
constant C{S) that only depends on S such that 



f^'ds\w"\^< r'ds((5|«;(3)|2 + c'((5)|w;f ) . 



(69) 



Choosing e.g. 5 = where C is the constant on the right hand side of eq. (68), and inserting 
into that equation, 



rSi 
J So 



ds|u;(^'|2 < C / ds 



So 



\wf <oo. 



(70) 



where C is independent of . Thus we may take the limit 5*1 — >■ oo in this equation, which shows 
y;(3) g L'^{So, oo). As a consequence of eq. (69), it follows that w" G L?{So, oo). Now we are ready 
to prove eq. (63) - we are going to show w, g{x{-), •), h{x{-), ■) e L^{Sq, oo). Firstly, for s > Sq, 



/•oo / /-oo J \ 1/2 / /.oo \ 1/2 

\w{s)\ < J \w'\ds < U ^] ij sV'dsj < Cs 



-3/2 
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where we have used Proposition 3 in the last inequahty. This estimate impHes w G L^(Sq,oo). 
This fact together with the following applications of Holder's inequality imply h{x, s) S L^{So, oo): 




Also, from eq. (58), «;W g L'^{So,(X}) for i e {0,1,2,3}, and g L^So,oo) nL°°(5o,oo) 
for i € {0,1}, it follows that g{x,s) G i-'^(S'o, oo), which in turn proves eq. (63). This con- 
cludes the proof of the proposition, since vA^^ = o(exp(— as)) as s — >■ oo for i = 0, 1,2,3 implies 
u{r),u'{r),'w{r),'w'{r) = o(exp(— ay^r))) as r — >■ oo. □ 

Proof of Theorem 1. First let A = 1. By Corollary 2, the functional is well defined and bounded 

from below. Let {(%, Wj)}jeN C W bo a minimizing sequence for E. By Corollary 2, E(uj,Wj) = 
E(uj,Wj), and {(uj,Wj)}j^\s4 is a minimizing sequence for E. By Proposition 2, it has a limit 
(uq, wq) e W and E{uo, wq) = infyy E = inf^ E, which implies E{uq, wq) = infw E. By Proposi- 
tion 4, 

uo{r) =^ + o(cxp(-CTVf)) 
Wo{r) =1 -f o(exp(-(TV^)) 

for all < (7 < 2^/2 as r — > oo. This proves the claim for A = 1. For general A, the functional is still 
bounded from below and well defined by by the change of variables eq. (11) with R = oo. Consider 
aminimizing sequence {{uj,Wj)}j(z^ oiEx. Then {{X^^Uj{X-),Wj{\-)}j^tt^ is aminimizing sequence 
for E, again by the change of variables eq. (11), and all the claims of Theorem 1 follow. □ 
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